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SUMMARY 

In this study, Coxiella burnetii seroprevalence was assessed for dairy and non-dairy sheep farm 
residents in The Netherlands for 2009-2010. Risk factors for seropositivity were identified for 
non-dairy sheep farm residents. Participants completed farm-based and individual questionnaires. 
In addition, participants were tested for IgG and IgM C. burnetii antibodies using 
immunofluorescent assay. Risk factors were identified by univariate, multivariate logistic 
regression, and multivariate multilevel analyses. In dairy and non-dairy sheep farm residents, 
seroprevalence was 66-7% and 51-3%, respectively. Significant risk factors were cattle contact, 
high goat density near the farm, sheep supplied from two provinces, high frequency of refreshing 
stable bedding, farm started before 1990 and presence of the Blessumer breed. Most risk factors 
indicate current or past goat and cattle exposure, with limited factors involving sheep. Subtyping 
human, cattle, goat, and sheep C. burnetii strains might elucidate their role in the infection risk 
of sheep farm residents. 
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INTRODUCTION 

Q fever, caused by Coxiella burnetii, is a worldwide 
zoonosis with goats, sheep, and cattle as primary 
sources for human infections [1]. Humans are usually 
infected by inhalation of contaminated aerosols 
originating from parturient animals and their birth 
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products [1-3]. Acute Q fever presents itself as a self- 
limiting febrile illness, pneumonia or hepatitis, with 
a small proportion developing chronic infections 
(mainly endocarditis and vascular infections) [4, 5]. 

From 2007 until 2009, large Q fever outbreaks 
occurred in The Netherlands, with over 3500 human 
cases notified [6]. Abortion waves at dairy goat 
farms were the primary source of these infections 
[7-9]. Between 2006 and 2008, C. burnetii abortion 
waves occurred on two dairy sheep farms [9]. 
Infected non-dairy sheep farms were not associated 
with an increased number of human cases living 
near these farms [10], although cases occurred in 
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individuals living a small distance from or having 
direct contact with non-dairy sheep in The 
Netherlands [11, 12]. Internationally, several sheep- 
related Q fever outbreaks have been reported 
[13-19]. 

In The Netherlands, sheep farms can be distin- 
guished from dairy farms and fat lamb-producing 
farms. There is a small dairy sheep industry with 
<50 farms, in which sheep are usually milked twice 
a day during several months each year. The number 
of sheep per farm differs from <50 to almost 1000 
with most kept outdoors for part of the year. On the 
fat lamb-producing sheep farms the sheep are kept 
outside, except for a few weeks around lambing, 
which usually occurs inside. Except for meat pro- 
duction, non-dairy sheep are also kept for breeding 
purposes or nature management. 

So far, no international studies have addressed 
the seroprevalence and risk factors for acquisition of 
C. burnetii infection in sheep farmers and their house- 
hold members. Therefore, our aim was to determine 
the C. burnetii seroprevalence in both dairy and non- 
dairy sheep farmers and their household members, 
and for the large non-dairy sector, to identify individ- 
ual and farm-related risk factors for seropositivity. 

MATERIAL AND METHODS 

All dairy sheep and non-dairy sheep farms in The 
Netherlands with at least 100 breeding ewes in 
November 2008, according to the national identifi- 
cation and registration database, were ehgible. A 
minimum of 100 ewes, considered to be a professional 
farm, was chosen because in the early stage of the 
Dutch epidemic it was clear that only (relatively 
large) commercial (dairy goat) farms were incrimi- 
nated as a potential source; no obvious role for 
small farms was observed [9]. Besides, smaller hobby 
farms have different management and farm residents 
of those farms are assumed to have a more limited 
exposure to sheep-related pathogens compared to 
commercial farms. Between September and 
December 2009, 32 dairy sheep farmers were ap- 
proached for the study. In addition, in March and 
April 2010, 1344 non-dairy sheep faimers were 
approached for participation. At the time of inclusion 
in 2010, those farms with at least 60 unvaccinated 
breeding animals were kept in the study. Farms 
with vaccinated sheep were excluded because in 
this integrated human-veterinary study the sheep at 
these farms were likely to be seropositive due to 



vaccination; vaccine-induced and naturally induced 
seroresponses cannot be distinguished to assess the 
true seroprevalence from natural infection. Second, 
we assumed that the infection rate for farm residents 
could be different for farms with vaccinated sheep 
(leading to reduced exposure) compared to farms 
with unvaccinated sheep. About 3 weeks after the 
initial invitation, all non-responding farmers were 
sent a written reminder. Because of the small number, 
dairy sheep farmers who did not respond to this 
second invitation were contacted by telephone. 

After written informed consent, a maximum of 
three persons were selected from each farm, i.e. the 
farmer and a maximum of two family members 
aged ^ 12 years residing at the farm; in some instances 
other persons working or living on the farm were 
selected. Each participant received a questionnaire 
addressing individual-based risk factors like age, 
gender, profession, ownership or contact with rumi- 
nants and pets, consumption of unpasteurized milk, 
medical history, and contact with agricultural pro- 
ducts. In addition, the farm owner or farm manager 
completed a farm-based questionnaire addressing 
characteristics like farm hygiene and management, 
herd size, presence of other hvestock and pets, stable 
environment, and lambing season characteristics. 
Separate farm-based questionnaires were developed 
for dairy farms and non-dairy farms because of clear 
differences in farm management. A professional lab- 
oratory assistant visited the farms to collect blood 
samples from all participating individuals for sero- 
logy. All data of the dairy sheep farms were collected 
between September 2009 and September 2010, for the 
non-dairy sheep farms data were collected between 
April and September 2010. The Medical Ethical 
Commission of the University Medical Center 
Utrecht approved the study protocol (no. 09-1 89/K). 

Serological analysis 

Serum samples were tested for C. burnetii IgM and 
IgG antibodies, both phases I and II, using an indirect 
immunofluorescence assay (IFA) with a screening 
dilution of 1:32. Participants without any positive 
antibody result and participants with a solitary IgM 
phase I or solitary IgM phase II result were classified 
as seronegative. All other outcomes were classified 
as seropositive. Those with IgM phase II antibodies 
were designated as 'relatively recent infections' and 
included possible current infections. The term 'rela- 
tively recent' was chosen as IgM phase II is found 
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to persist in the majority of cases for 1 year 
post-infection and may even persist up to 4 years 
post-infection [20, 21] (C. C. H. Wielders, personal 
communication). Seropositives without IgM phase II 
antibodies were designated as 'past infections'. As 
the latter group also includes possible chronic infec- 
tions, a further distinction was made between serologi- 
cal profiles that had IgG phase I ^ 1:1024 indicative 
for a chronic infection according to the new Dutch 
consensus guidelines [22]. 

Statistical analyses 

Dairy sheep farms 

All data were analysed with SAS, version 9.2 (SAS 
Institute Inc., USA). For the dairy sheep farms in 
The Netherlands, participation bias was investigated 
by comparing participating and non-participating 
farms with regard to herd size, urbanization degree 
and region. The seroprevalence of C. burnetii in resi- 
dents and the corresponding 95% confidence interval 
(CI) were calculated. Descriptive statistics were per- 
formed by analysing frequency tables and studying 
distributions of continuous variables. No risk factor 
analysis was performed because of the small number 
of participants. 

Non-dairy sheep farms 

To study participation bias, participating and non- 
participating farms were compared with regard to 
herd size, cattle, sheep, and goat density in the sur- 
roundings, urbanization degree, region, situated inside 
or outside a compulsory Q fever vaccination area, 
number of bulk-milk-positive dairy goat or dairy 
sheep farms in a radius of 5 and 10 km, and distance 
in metres to the closest bulk-milk-positive small 
ruminant farm. 

The seroprevalence of C. burnetii and the corre- 
sponding 95% CI were calculated. For descriptive 
statistics, frequency tables were analysed. In addition, 
distributions of continuous variables were studied, 
and if not linearly related to the outcome variable, 
continuous variables were recoded into classes. 

Univariate logistic regression analysis was per- 
formed to assess the main factors associated 
with C. burnetii seropositivity at the individual level 
[P<0-20 in the hkelihood ratio test (-2LL)]. 
Variables with <20 participants in one risk category 
were excluded. Age was always kept in the model 
because of the frequent association with Q fever 



seropositivity in the literature. Proxy outcomes, such 
as sheep seropositivity, were not included in the multi- 
variate analysis. If several variables, which were 
associated in the univariate analysis, were interrelated, 
a preferred variable was chosen and related variables 
were excluded. The preferred variable was chosen 
based on the most informative value, the strongest 
association or most relevant exposure (exposure at 
own farm instead of comparable exposure at other 
farms). All identified individual variables were ana- 
lysed with a manual backwards elimination procedure 
until all variables were significant at the 10% signifi- 
cance level in the likelihood ratio test, starting with 
a full multivariate logistic regression model. 

Subsequently, potential risk factors derived from 
the farm-based questionnaire were analysed by uni- 
variate multilevel analyses considering clustered farm- 
based data for all persons within the same farm, using 
a unique fami number as cluster variable. All farm 
variables which were significant in the univariate 
analysis (/'<0-20), were analysed with a manual back- 
ward elimination procedure starting with a full multi- 
level model. 

Finally, both the individual and farm-based charac- 
teristics from the two final submodels were combined 
in a multivariate multilevel analysis to identify the 
independent risk determinants for C. burnetii sero- 
positivity. The final model fit was assessed by the 
quasi-likelihood under the independence model 
criterion (QIC) goodness-of-fit statistic for generalized 
estimation equation (GEE) models. 



RESULTS 
Dairy sheep farms 

Out of the 32 invited farms, 12 participated (response 
rate 37-5%). The participating farms were all situated 
in a rural area (<500 addresses/km^). Participating 
and non-participating farms were comparable with 
regard to urbanization degree and province distri- 
bution. However, participating farms had a median 
number of 529 sheep (range 143-1163) v.v. the signifi- 
cantly lower median of 353 sheep (range 96-730) for 
non-participating farms (/'=0-03). 

Twenty-seven study participants (mean age 
38-7 years, range 14-61, 63% male), provided a 
blood sample. Overall, 18 (66-7%) participants were 
seropositive: 80-0% for the 15 farmers (12 males), 
and 50-0% for the 12 household members (five 
children, five female spouses, one male spouse, one 
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seasonal worker). Three (11-1%) participants had a 
relatively recent C. burnetii infection (IgM phase II 
antibodies). None consulted their general practitioner 
or were hospitalized because of influenza-like illness 
or fever. One participant had an IgG phase I titre of 
^ 1:1024, indicating a possible chr^onic case [22]. 



Non-dairy sheep farms 

Non-response analyses 

Out of the 1344 approached farms, at least 32 
appeared to be no longer eligible because they had 
<60 animals at inclusion or had vaccinated all their 
sheep. Of the remaining 1312 farms, 119 participated 
in the study (response rate 9T%). 

A significant difference was found for sheep 
density in the 5-km radius of participating and non- 
participating farms, 34-5 (range 1 ■8-143-6) and 47-5 
(range 1-0-162-9) sheep/km'^ in the 5-km radius 
(excluding own sheep), respectively (f = 0-01). In 
addition, the number of sheep was borderline signifi- 
cantly higher at the participating farms (median 
191 sheep, range 102-1310), compared to the non- 
participating farms (median 167 sheep, range 
100-2857). For the other variables, no significant 
differences were found between participating and non- 
participating farms (Table 1). 

Descriptive characteristics 

The 119 participating farms were mainly situated 
in the provinces of Noord-Holland and Friesland, 
commonly (90-8%) situated in rural areas (<500 
addresses/km^) and the most common breeds at 
the farms were Texel (57-0%) and Swifter (46-5%). 
The farms were mainly started after 1950 (9-6% 
1875-1950, 39-4% 1951-1980, 51-0% after 1980). 
Out of the 1 14 farms with a farm-based questionnaire, 
23 (20-2%) kept one or more goats, 45 (39-5%) kept 
dairy cattle and/or beef cattle, and 13 (11-4%) other 
farms reported that cattle were present on their pas- 
tures. The farms could have one or more function; 
95 (83-3%) farms kept sheep for meat production, 
53 (46-5%) farms for rearing, and 20 (17-5%) farms 
for nature management. Of those 20 farms, 12 farms 
kept their sheep exclusively for nature management. 

From the 119 farms, 271 persons provided a blood 
sample (mean age 47, range 12-93 years, 55% male). 
Of those, 266 completed the individual self- 
administered questionnaire and from 261 individuals 



information was available from the farm-based 
questionnaire. 

C. burnetii seroprevalence was 51-3% (95% CI 
45-5-57-4). In the univariate analysis, seroprevalence 
was significantly higher for farmers (58-8% vs. 36-3% 
for spouses) and for males (57-7% vs. 43-4% for 
females). Out of the 271 participants, seven (2-6%) 
had a relatively recent infection (IgM phase II anti- 
bodies). No participant had an IgG phase I titre sug- 
gestive for chronic infection. 

Although the seroprevalence of the farm residents 
was higher for those hving on a dairy sheep farm, 
the difference was not statistically significant [odds 
ratio (OR) 1-9, 95% CI 0-8^-4] for dairy sheep farm- 
ers vs. non-dairy sheep farmers). 

Univariate analyses at individual and farm level 

All individual and farm-based variables, which were 
tested in the univariate analysis for relationship with 
human C. burnetii seropositivity, are displayed in 
Tables 2 and 3. 

Multivariate and multilevel analyses 

In the multivariate analyses, from 23 individual vari- 
ables which were associated in the univariate analysis, 
four were independently associated with C. burnetii 
seropositivity (Table 4). In addition, 10/23 farm-based 
variables included in the multilevel analyses were sig- 
nificantly independent risk or protective factors and 
together were used as the full multilevel start model 
(Table 5). 

Combined multilevel analyses of individual and 
farm-based factors 

In the final combined multilevel model, significant risk 
factors were contact with cattle at own or other farm, 
past employment in the cattle sector, high goat density 
in the vicinity of the farm, living or working at a 
farm that was started in 1990 or later, the presence 
of Blessumer breed on the farm, cattle on the same 
pastures used by sheep, although not simultaneously 
with the sheep, high frequency of refreshing the bed- 
ding in the sheep stables, and sheep supplied from 
the provinces of Groningen or Noord-Holland 
(Table 6). Borderline significant risk factors were age 
40^9 years, and presence of dairy cattle during the 
stabling period of the sheep. In addition, sheep lamb- 
ing outside was a significant protective factor, and air 
entering the stable through the door was a borderline 
significant protective factor. 
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Table 1 . Non-response analyses of non-dairy sheep farms, comparison of participating and non-participating farms 



Numerical variables 


Participating 
farms (A' =119) 
Median 




Non-participating 
farms (iV=1193) 
Median 


P value 


Number or sheep 
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Sheep density (number of sheep/km^ excluding own animals 


34-5 




47-5 


0-01 


in a 5-km radius)* 










Closest Q fever bulk-milk-positive dairy goat or dairy 


13960 




13 806 


0-70 


sheep farm (metres)* 










Number Q fever bulk-milk-positive dairy goat or dairy 


0 (min = 0, max = 
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Categorical variables 


n (%) 




M (%) 


P value 


Inside vaccination area 


20 (16-8) 




181 (15-2) 


0-64 


Outside vaccination area 


99 (83-2) 




1012 (84-8) 




Urbanization 










Very high urban area*t 


0 (0 0) 




2 (0-2) 


0-37 


Higli urban area 


0 (0 0) 




3 (0-3) 




Moderate urban area 


4 (3-3) 




14 (1-2) 




Minor urban area 


7 (5-9) 




84 (7-0) 




Rural area 


108 (90-8) 




1086 (91-3) 




Province 










Drenthe* 


4 {i 4) 




3/ [^0) 


U 3Z 


r icvuiaiiii 


1 KVi 6) 
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14/11 .B\ 
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1 70 n zl-'^^ 
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Groningen 


11 (9-2) 




93 (7-8) 




Limburg 


4 (3-4) 




23 (1-9) 




Noord-Brabant 


12 (101) 




74 (6-2) 




Noord-Holland 


29 (24-4) 




241 (20-3) 




Overijsel 


11 (9-2) 




86 (7-2) 




Utrecht 


2(1-7) 




48 (4-1) 




Zeeland 


2(1-7) 




49 (4-1) 




Zuid-Holland 


11 (9-2) 




126 (10-6) 





A', Total number of individuals. 

* Four missing values at non-participating farms. 

t Urbanization degree: very high urban area >2500 addresses/km^; high urban area = 1500-2500 addresses/km"; 
moderate urban area= 1000-1500 addresses/km^; minor urban area = 500-1000 addresses/km^; rural area <500 addresses/km". 



community in the epicentre of the Q fever outbreak 
in 2007 (25-1%), and in blood donors in the most Q 
fever-affected areas in 2009 (12-2%), indicating that 
sheep farm residents have an increased life-time risk 
of acquiring a C. burnetii infection compared to the 
general Dutch population [7, 23, 24]. 

The observed seroprevalence in Dutch sheep 
farm households is also high compared to a study 
of sheep farmers in Sweden (28-5%) [25], and of 



DISCUSSION 
Seroprevalence 

The seroprevalence of non-dairy (51-3%) and dairy 
sheep farm residents (66-7%) is clearly higher com- 
pared to the seroprevalence estimate of 2-4% in the 
general population before the outbreak occurred 
in The Netherlands in 2006-2007. It is even higher 
compared to the seroprevalence found in a small 
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Table 2. Univariate logistic model of individual factors related to C. burnetii seropositivity in non-dairy sheep farm 
residents (P<0-20, -2LL) 







Frequency 


Sero- 








(AO 


prevalence 




Variable 


Category 


{N=26€) 


(%) 


OR (95% CI) 


Gender* 


Male 


144 


57-6 


l-77(l-09-2-88) 




Female 


122 


43-4 


Reference 


Age (years)* 


12-19 


21 


57-1 


2-04 (0-72-5-76) 




20-39 


45 


511 


l-60(0-70-3-63) 




40-49 


68 


58-8 


2-18 (1-03-4-63) 




50-59 


84 


50-0 


1-53 (0-74-3-13) 




>60 


48 


39-6 


Reference 


Work and/or live on farm 


Work and live 


188 


53-7 


1-61 (0-83-3-15) 




Work, but not live 


35 


48-6 


1-31 (0-53-3-22) 




Not work, but hve 


43 


41-9 


Reference 


Function 


Farmer 


136 


58-8 


2-51 (1 •42^-44) 




Spouse 


80 


36-3 


Reference 




Childf 


39 


53-9 


2-05 (0-94-4-46) 




Otherl 


11 


54-6 


2-11 (0-59-7-53) 


How often in stable 


Every day 


185 


55-7 


Reference 




Every week 


56 


41-1 


0-56 (0-30-1 -02) 




Every month 


10 


500 


0-80 (0-22-2-84) 




Less than once a month/never 


15 


33-3 


0-40 (0-13-1-21) 


Amount of work at farm* 


Full working week 


61 


63-9 


2-39(1 ■25^-56) 




Up to half a working week 


97 


52-9 


1-49(0-86-2-59) 




Never/occasionally 


108 


42-6 


Reference 


Feed sheep* 


Yes 


225 


55-6 


3-41 (1-63-7-14) 




No 


41 


26-8 


Reference 


Load and imload sheep 


Yes 


194 


56-2 


2-14 (1-23-3-72) 




No 


72 


37-5 


Reference 


General healthcare of sheep 


Yes 


201 


55-7 


2-15(1-21-3-82) 




No 


65 


36-9 


Reference 


Remove manure 


Yes 


180 


57-8 


2-31 (1-36-3-92) 




No 


86 


37-2 


Reference 


Spread manure* 


Yes 


124 


58-9 


1-80 (1-10-2-92) 




No 


142 


44-4 


Reference 


Clean stables 


Yes 


167 


56-3 


1-75(1-06-2-89) 




No 


99 


42-4 


Reference 


Administrative work 


Yes 


193 


54-4 


1-62 (0-94-2-78) 




No 


73 


42-5 


Reference 


Wear overalls or boots* 


Yes 


234 


54-3 


3-03 (1-35-6-84) 




No 


32 


28-1 


Reference 


Contact with cattle at own or 


Yes 


172 


63-4 


4-29 (2-49-7-40) 


other farm*§ 


No 


94 


28-7 


Reference 


Contact with horses at own or 


Yes 


145 


59-3 


2-07(1-27-3-38) 


other farm*§ 


No 


121 


41-3 


Reference 


Contact with pigs at own farm*§ 


Yes 


24 


37-5 


0-54(0-23-1-29) 




No 


242 


52-5 


Reference 


Indirect contact with poultry 


Yes 


93 


57-0 


1-44(0-87-2-39) 


at own farm* II 


No 


173 


48-0 


Reference 


Indirect contact with rats 


Yes 


45 


64-4 


1-93 (0-99-3-76) 


at own farm* II 


No 


221 


48-4 


Reference 


Contact with goats at other farm*§ 


Yes 


32 


62-5 


1-70 (0-79-3-63) 




No 


234 


49-6 


Reference 


Contact with sheep at other farm*§ 


Yes 


102 


60-8 


1-89(1-14-3-12) 




No 


164 


45-1 


Reference 


Contact with dogs at other farm*§ 


Yes 


112 


58-9 


1-72(1-05-2-82) 




No 


154 


45-5 


Reference 
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Table 2 (cont.) 







Frequency 


Sero- 








(AO 


prevalence 




Variable 


Category 


(Af=266) 


(%) 


OR (95% CI) 


Indirect contact with poultry 
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1 ou (u y J— z I z) 
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Direct contact with wool* 
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Li J 


ou z 
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Reference 


Direct contact with hay, straw 


Viae 

I es 


HQ 


34 o 


9.QQ (\ .A 1 A.9QA 

z 70 (^1 41— 0 zy^ 


or animal feed* 


rvO 


1Q 

JO 


zv u 


Reference 


Direct contact with raw milk 


Viae 
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1 7 1 \ V lU— 3 3Z ) 
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Vyj 
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Reference 
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uu / 


9-1 7 n -1 1—4. 9A^ 
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Oil 

ZZV 


4o U 


Reference 


Direct contact with cattle manure 


I es 


1 1 n 


Do Z 


T.TA (\ .07 C.C7\ 

3 3U (1 7/— J jz^ 
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i J J 




Reference 


Direct contact with live-born animals 


I es 






"X-Al (■! .9 1 Q.AQ\ 
3 ^Z \ V Zl— 7 07^ 


UUllllg lalllUlIlg pcilUCl 


IN U 


zu 


ZJ u 
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Direct contact with dead-born 


I es 


9 1 n 
ziu 


j4 J 


1 o4 (^1 Ul— 3 3Dj 


animais/ pidcen id 


XTn 


jd 


37 J 


IvCiCICllCC 


Xick bite* 


Viae 
I cS 


01 


'^-Z 0 


U 0^- (^U 30—1 V^) 




ISO 


zUj 


JJ / 


Reference 


Did not work in animal 


Viae 

Yes 


1 1 A 
i l-^ 


3y J 


U 44 \\j Z i—\J 1 Z) 


husbandry/agriculture in the past 


iNO 


VjZ 




Reference 


Employment in cattle sector in the past* 


I es 


1 n7 
iU/ 


04 J 


Z 47 (1 jU— 4 14J 




No 


159 


42-1 


Reference 


Worked in animal transport/ transport 


Yes 


37 


70-3 


2-56(l-21-5-42) 


of agricultural products in the past* 


No 


229 


48-0 


Reference 


As a child lived at: 


Cattle farm 


151 


59-6 


2-04(M8-3-53) 




Other kind of farm 


34 


35-3 


0-75 (0-33-1-73) 




No farm 


81 


42-0 


Reference 


As a child worked in animal care/with 


Yes 


178 


56-2 


l-85(M0-3-ll) 


manure/hay/in vegetation care* 


No 


88 


40-9 


Reference 



A', Total number of individuals; OR, odds ratio; CI, confidence interval, — 2LL, likelihood ratio test. 
* Variables included in subsequent multivariate individual analyses before manual backward elimination. 

I Children aged <18 years («= 17) and older children (« = 22) of the farmer. 
{ Employees, shepherds, other family members. 

§ See animals at <5 m or touch animals. 

II See animals at <5 m. 



farmers from all types of farms: 17-8% in Poland, and 
27-3% in the UK [26, 27]. Generally, it is difficult 
to compare international seroprevalence studies, 
because most studies use different tests or cut-off 
values. The cut-off value of the test in our study 
(^1:32) was chosen because it allowed comparison 
with other population surveys conducted in The 
Netherlands [23, 28]. 

Dairy sheep farm residents had a higher seropreva- 
lence compared to non-dairy sheep farm residents. 
Although no statistically significant difference in sero- 
prevalence was found between the residents of both 



farm types, this might be due to lack of power because 
of the small number of participants from dairy 
sheep farms. In this study it was impossible to assess 
which risk factors were responsible for the higher sero- 
prevalence in dairy sheep farm residents, due to the 
low number of participating dairy sheep farm resi- 
dents. In addition, because of the differences in farm 
management, the farm-based questionnaires of both 
farm types were not the same, therefore pooling 
the analysis with the other sheep farm residents 
to increase power was not an option. Specific 
research, targeting all current dairy sheep farms 
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in The Netherlands (7i~40), might elucidate further 
risk factors next to the higher sheep seroprevalence, 
explaining the higher seroprevalence in dairy sheep 
farm residents. Nevertheless, it might well be that 
dairy farm residents were more exposed to Coxiella, 
as the seroprevalence in dairy sheep at these same 
farms was signiticantly higher compared to that of 
non-dairy sheep (data not shown). A higher vulner- 
ability for infection of breeds selected for milk 
production rather than for disease resistance has pre- 
viously been observed for dairy cattle, dairy sheep, 
and dairy goats [29, 30]. In addition, dairy sheep 
are more often housed in stables compared to non- 
dairy sheep which spend most of the year outside. 
Indoor housing might facilitate the spread of C. bur- 
netii in dairy sheep and to humans. Moreover, the 
higher seroprevalence in dairy farm residents 
might be explained by more intense contact with 
dairy sheep. 

The seroprevalence of the dairy sheep farm resi- 
dents (66-7%) was comparable to the seroprevalence 
of dairy goat farm residents (68-7%) in The 
Netherlands [28]. Furthermore, the percentage of rela- 
tively recent infections (clinical status unknown as no 
questions addressed current Q fever compatible symp- 
toms) in the dairy sheep farm residents (11-1%) is 
comparable to that of the dairy goat farm residents 
(11-2%) [28]. Additionally, the percentage of partici- 
pants with an indication for a possible chronic infec- 
tion is also similarly high for dairy sheep and dairy 
goat farm residents (3-7% and 4T%, respectively) 
[28]. In contrast, the percentage of relatively recent 
infections and possible chronic infections are lower 
for non-dairy sheep farm residents (2-6% and 0%, 
respectively). Therefore, currently C. burnetii infection 
seems to be a more serious and on-going health prob- 
lem in dairy goat and dairy sheep farm residents 
compared to non-dairy sheep farm residents, although 
the numbers are relatively small. 

Although numbers are too low to draw any con- 
clusion and do not allow for valid statistical testing, 
the 10 (three from dairy and seven from non-dairy 
farms) relatively recent (IgM phase II positive) cases 
were generally younger (median 37 years v.y. median 
50 years for the seronegatives), were more often 
male (80% vs. 48%) and more often lived on a dairy 
sheep farm (30% of the recently infected vs. 6% 
of the seronegatives). This may point to ongoing 
infections especially in male dairy sheep farm re- 
sidents, in the relatively early days of their contact 
with sheep. 



Risk and protective factors for non-dairy sheep 
farm residents 

One of the protective factors for C. burnetii seroposi- 
tivty was sheep lambing outside. Farm residents 
might be less exposed to contaminated aerosols in 
that situation, compared to lambing inside stables. 

In addition, several risk factors for C. burnetii sero- 
positivity were identified in this study. McCaughey 
et al. [31] suggested in his study in the general popu- 
lation (age 12-64 years) that most people acquired 
C. burnetii infection between ages 25 and 34 years 
and after that age seroprevalence remained stable. 
This age trend was not seen in our study; sheep farm 
residents had already a high seroprevalence at young 
age (12-19 years). This might be explained by 
exposure to infected animals at a young age. The high- 
est seroprevalence found in humans (age 40^9 years), 
matches the most common age group of notified clini- 
cal Q fever cases in The Netherlands [9]. The increased 
risk at this age seems not to be explained by differ- 
ences in specific work activities, frequency of cattle 
contact, or hours worked. Perhaps host factors play 
a role in the increased risk, or it generally reflects reg- 
ular exposure to the bacterium and repeated develop- 
ment of antibodies (booster effect), not adequately 
measured by the questions in the questionnaire. 

Animal movement is a known risk factor for the 
transfer of microorganisms and should be discour- 
aged [32, 33]. Why specifically supply of sheep 
from the northern provinces of Noord-Holland 
and Groningen showed an independent increased 
risk for infection of the farm residents is not 
clear. The seroprevalence in sheep in these two pro- 
vinces was not significantly different from pre- 
valences in other provinces, both in the current 
study (B. Schimmer et al., unpublished data) and in 
a previous study in 2008 using convenience serum 
samples from sheep [30]. 

It is also unknown why the fact that a farm started 
before 1990 was a risk factor. No change in farm man- 
agement is known around that year that could influ- 
ence the risk of a C. burnetii infection. 

Having the Blessumer sheep breed on the farm 
was the next significant risk factor. This breed is a 
crossing of the breeds of Texel (non-dairy sheep) 
and Flemish sheep (dairy sheep); therefore, the 
Blessummer breed might have a lower disease resist- 
ance [29, 30]. Differences in infection rates between 
sheep breeds have not yet been studied to investigate 
whether Blessumer sheep are more often infected. 
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Table 3. Univariate multilevel analysis of farm-based factors related to C. burnetii seropositivity in non-dairy sheep 
farm residents (P<0-20) 







Frequency 


Sero- 








(AO 


prevalence 




Variable 


Category 


(7V=261)* 


(%) 


OR (95% CI) 


Urbanizationt{§ 


Moderate or minor 
urban area 


28 


67-9 


2-00 (0-80-5-04) 




Rural area 


242 


49-2 


Reference 


Goat density (number of goats/km excluding 


<2-9 


135 


38-5 


Reference 


own animals in a 5-km radius)t§ 


2-9-1 1-3 


67 


68-7 


3-59 (1-86-6-91) 




^11-4 


68 


58-8 


2-38 (1-18^-79) 


Sheep density (mmiber of sheep/km excludmg 


<33-7 


133 


41-4 


Reference 


own animals in 5-km radius)t§ 


33-7-79-0 


69 


53-6 


1-68 (0-87-3-25) 




^79-1 


68 


67-7 


2-98 (1-54-5-78) 


Cattle density (number of cattle/km in the 


< 200-0 


240 


47-9 


Reference 


municipality) t§ 


^200-0 


30 


76-7 


3-20 (1-37-7-51) 


Number of Q fever bulk-milk-positive dairy 


0 


166 


45-8 


Reference 


goat or dairy sheep farms in a 10-km radiust§ 


1^ 


104 


59-6 


1-78 (1-02-3-11) 


Closest Q fever bulk-milk-positive dairy goat 


<5-0 


35 


62-9 


Reference 


or dairy sheep farms (km)§ 


5-0-9-9 


69 


58-0 


0-39 (0-14-1-13) 




10-0-14-9 


53 


41-5 


0-87 (0-30-2-54) 




15-0-19-9 


41 


61-0 


0-82 (0-32-2-14) 




^ 20-0 


72 


40-3 


0-42 (0-16-1-10) 


Year farm startedf 


Before 1990 


165 


44-2 


Reference 




1990 or later 


75 


61-3 


1-97 (1-12-3-48) 


Distance between house and pastures 


<30m 


127 


40-2 


Reference 




^30m 


103 


61-1 


2-20 (1-23-3-94) 


Nimiber of male sheep 201 Of 


<6 


130 


46-9 


Reference 




6-20 


56 


60-7 


1-78 (0-85-3-75) 




>20 


41 


51-2 


1-20 (0-53-2-70) 




No 


16 


56-3 


1-30 (0-42^-00) 


Zwartbles breed present on farmf 


Yes 


30 


63-3 


1-75 (0-89-3-42) 




No 


228 


48-7 


Reference 


Rijnlam breed present on farm 


Yes 


7 


85-7 


5-72 (0-78^2-12) 




No 


251 


49-4 


Reference 


Blessumer breed present on farmf 


Yes 


21 


76-2 


3-51 (1-25-9-81) 




No 


237 


48-1 


Reference 


Animals at same pasture simultaneously 


None 


160 


52-5 


Reference 


with sheep 


Cattle 


66 


59-1 


1-30 (0-73-2-33) 




Other 


27 


18-5 


0-21 (0-07-0-66) 


Cattle at same pasture but not 


Yes 


62 


74-2 


3-90 (1-74-8-72) 


simultaneously with sheepf 


No 


188 


42-0 


Reference 


Straw bedding in the stables 


Yes 


243 


50-2 


0-69 (0-40-1-21) 




No 


5 


60-0 


Reference 




No stable 


10 


50-0 


0-31 (0-24-1-68) 


How often bedding in stable 


Every other day or more 


200 


53-0 


1-77 (0-83-3-76) 


is refreshedf 


Once or tvsace a week 


47 


38-3 


Reference 




No stable 


10 


50-0 


1-46 (0-49^-35) 


Air enters stable through doorf 


Yes 


163 


46-6 


0-64 (0-35-1-18) 




No 


79 


58-2 


Reference 




No stable 


10 


50-0 


0-67 (0-25-1-80) 


No farm animals present on farm 


Yes 


73 


42-5 


0-63 (0-34-1-14) 


other than sheep 


No 


183 


53-6 


Reference 


Other farm animals present in 


Yes 


164 


54-9 


1-71 (0-98-3-00) 


sheep stables 


No 


92 


42-4 


Reference 


Laying hen in stablef 


Yes 


35 


65-7 


2-11 (0-88-5-04) 




No 


215 


47-9 


Reference 
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Table 3 (cont.) 







Frequency 


Sero- 








(AO 


prevalence 




Variable 


Category 


(A^=261)* 


(%) 


OR (95% CI) 


Dairy cattle in stabley 


1 es 


00 


71 .9 
/ 1 Z 


DDI \ V 1 D— D 43 ) 




INO 


1 9.A 


4Z y 


Reference 


Type ol leed method 


By hand/ wheelbarrow 


one 

ZUo 


/1ft. 1 


Reference 




iviixer 


1 /I 
14 


7 1 .A 
/ 1 4 


z yi yz— y z 




Shovel 


'X'X 


/IS.'n 
4o J 


1 .09 i^n.^l 1 .Q7A 
1 uz \\j 3 J— 1 y 1 ) 


l-,anibing outsidey 


I es 


11 

Z / 


J 1 \) 


U 33 \\J ZD— 1 ZAj) 




INO 




Z>l J 


Reference 


Number ol yearlings whicn lambed 




zUcS 


40 0 


Reference 


m zuuvT 


/in 




oz u 


1 .70 i^n.SQ T./^TA 
V ly yJ oy—D KiD) 


Number dead-born lambs in 2009 




4y 


4U o 


Reference 




0—14 




J / u 


1 .sc /^n.c^ /1 . 1 ^\ 

1 oo (^U o3— 4 13^ 




1 J — Z'f 




410 


1 .no ((\.Ai 9.^n\ 
1 uy (U 4/— z 3Uj 




'-'Zj 


fo 


J4 Z 


1 \jy yJ 1 1 — 4 xjj) 


/VDoriion rdie zuu/, zuuo, zwuy^^/ojj 


*-^4 in all three years 


1 

ly J 


40 Z 


Reference 




^4 in at least one year 


C 1 

j1 


DO* / 


7.Q ^ /'I . 1 7 /1 .07\ 
ADD (1 IZ— 4 yZj 


Afterbirth of normally lambed animalf 


Leave in stable or pasture 


50 


580 


Reference 






o4 




U 04 \\J D\J—l 30^ 




Direct or once a day 


1 nn 

lUU 


jl U 


n.79 ((\.'i.A 1 .^"^A 
U /Z (^U 34— 1 33 j 




render bucket 










111 ■f^^^^'f rw rw\ r^p^ o 

uiiccL ui uiicc d. uay 


zu 




u 3 i yj 1 u — u y 1 ) 




manure yard 


















rarm tenure j 


Closed for ewes and rams 

Ul Uiliy LIUSCU iUl CWCa 


1 


43 'Z 


Reierence 




Not closed for ewes 


79 

/z 


Oj J 


9.17 (^ .9/1 A.'\A\ 
L Dl ^1 Z4— 4 34^ 




allU lalllo 








Sheep supplied irom Oroningenj 


I es 


9A 

zO 


oy z 


Z 3U (U Sz— / 3 / ) 




INO 


ZZO 


/1Q.7 


Reference 


Sheep supplied Irom Noord- Brabantj 


I es 


97 

Z/ 


Oj\J 


1 .03 ((\.Si 1 /I . ^QA 
1 y3 (^U ol— 4 3oj 




InO 


99^ 

ZZ3 


4o y 


Reference 


Sheep supplied from Noord- Holland! 


Yes 


76 


59-2 


1-67 (0-89-3-15) 




No 


176 


46-6 


Reference 


Sheep supplied from Utrecht 


Yes 


15 


73-3 


2-69 (0-73-9-86) 




No 


237 


490 


Reference 


Presence of hygienic locker room 


Yes 


19 


68-4 


2-32 (0-81-6-62) 




No 


231 


48-5 


Reference 


Presence of disinfection bucketf 


Yes 


36 


611 


1-80 (0-89-3-65) 




No 


214 


48-1 


Reference 



A', Total number of individuals; OR, odds ratio; CI, confidence interval. 
* Not all numbers add up to the total due to missing values. 

t Variable included in later multivariate farm-based analyses before manual backward elimination. 

:|: Urbanization degree: moderate urban area = 1000-1500 addresses/km^; minor urban area = 500-1000 addresses/km^; 
rural area <500 addresses/km^. 

§ For the geographical data, information was available for all 270 individuals, including the nine people without a farm- 
based questionnaire. 



In the environment of dairy goat farms with a his- 
tory of abortion waves and of farms having 
PCR-positive bulk milk, relatively high levels of 
C. burnetii DNA were found [34]. A high goat density 
in the surrounding area of a participating farm is 
therefore considered a plausible risk factor for people 



living in the vicinity at the time of data collection. 
This was also demonstrated in several local outbreak 
investigations in The Netherlands in 2008-2009 [7, 8]. 

Maredly, several risk factors for C. burnetii sero- 
positivity in non-dairy sheep farm residents point to 
cattle exposure at present or in the past. This might 
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Table 4. Results of the multivariate logistic regression analysis for the individual characteristics (P<0-10, —ILL) 


in relation to non-dairy sheep farm residents C. burnetii seropositivity 




Variable 


Category 




Age (years) 


12-19 


2-81 (0-85-9-35) 




20-39 


1-42 (0-57-3-54) 




40-49 


2-29 (1-00-5-24) 




50-59 


1-12 (0-50-2-48) 




>60 


Reference 


Amount of work at farm 


Full working week 


2-42 (1-13-5-15) 




Up to half a working week 


1-23 (0-65-2-33) 




Never/occasionally 


Reference 


Contact with cattle at own or other farm* 


Yes 


3-87 (2-13-7-04) 




No 


Reference 


Worked in cattle sector in the past 


Yes 


1-79 (1-01-3-18) 




No 


Reference 


OR, Odds ratio; CI, confidence interval; — 2LL, likelihood ratio test; AIC, Akaike's Information Criterion. 


Number of observations used: 266 (AIC = 340-38). 






* See animals at <5 m or touch animals. 






Table 5. Results of the multilevel analysis with farm-based characteristics (P<0-10) as 


independent factors 


in relation to non-dairy sheep farm residents C. burnetii seroprevalence 




Variable 


Category 


OR (9S% CD 


Goat density (number of goats/km^ excluding own 


<2-9 


Reference 


animals in a 5-km radius) 


2-9-1 1-3 


1-60 (0-75-3-43) 




^11-4 


3-80 (1-67-8-65) 


Year farm started 


Before 1990 


Reference 




1990 or later 


3-97 (1-79-8-82) 


Blessumer breed present on farm 


Yes 


5-19 (2-36-11-41) 




No 


Reference 


Cattle at same pasture but not simultaneously with sheep 


Yes 


5-14 (2-17-12-19) 




No 


Reference 


How often bedding in stable is refreshed 


Every other day or more 


3-24 (1-49-7-07) 




Once or twice a week 


Reference 




No stable 


8-91 (2-17-36-68) 


Air enters stable through door 


Yes 


0-46 (0-23-0-92) 




No 


Reference 




No stable 


8-91 (2-17-36-68) 


Dairy cattle present during stabling period of sheep 


Yes 


3-33 (1-17-9-46) 




No 


Reference 


Lambing outside 


Yes 


0-34 (0-14-0-86) 




No 


Reference 


Sheep supplied from Groningen 


Yes 


4-17 (1-59-10-97) 




No 


Reference 


Sheep supplied from Noord-HoUand 


Yes 


3-93 (1-74-8-90) 




No 


Reference 



OR, Odds ratio; CI, confidence interval; QIC, quasi-likelihood under the independence model criterion. 
Number of observations used: 212. Number of levels used: 107 (QIC = 232-9560). 



suggest that cattle were partially responsible for the 
infections observed in the sheep farm residents. In a 
previous study in farmers (all farm types) contact 
with cattle was also described as a risk [27]. A recent 



published review including worldwide studies, 
suggested a higher seroprevalence of C. burnetii in 
cattle compared to goat and sheep [35]. In The 
Netherlands, a prevalence of 78-6% for antibodies 
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Table 6. Results of the multilevel analysis with individual and farm-based characteristics (P<0-10) as independent 
factors in relation to non-dairy sheep farm residents C. burnetii seroprevalence 



Variable 


Category 


OR (95% CI) 


Age (years) 


12-19 


0-96 (0-29-3-21) 




20-39 


1-96 (0-56-6-90) 




40^9 


2-43 (0-98-6-04) 




50-59 


1-54 (0-63-3-78) 




>60 


Reference 


Contact with cattle at own or other farm* 


Yes 


2-32 (1 •02-5-29) 




No 


Reference 


Worked in cattle sector in the past 


Yes 


3-98 (1-71-9-25) 




No 


Reference 


Goat density (number of goats/km excluding own 


<2-9 


Reference 


animals in a 5-km radius) 


2-9-1 1-3 


1-11 (0-46-2-68) 




^11-4 


5-86 (1-81-18-95) 


Year farm started 


Before 1990 


Reference 




1990 or Later 


3-67 (1-45-9-31) 


Blessumer breed present on farm 


Yes 


4-49 (1-59-12-65) 




No 


Reference 


Cattle at same pasture but not simultaneously with sheep 


Yes 


5-77 (2-29-14-56) 




No 


Reference 


How often bedding in stable is refreshed 


Every other day or more 


4-58 (1-69-12-37) 




Once or twice a week 


Reference 




No stable 


8-34 (1-71^0-60) 


Air enters stable through door 


Yes 


0-47 (0-21-1-01) 




No 


Reference 




No stable 


8-34(1-71^0-60) 


Dairy cattle present during stabling period of sheep 


Yes 


2-69 (0-81-8-95) 




No 


Reference 


Lambing outside 


Yes 


0-33 (0-12-0-92) 




No 


Reference 


Sheep supplied from Groningen 


Yes 


5-05 (1-73-14-69) 




No 


Reference 


Sheep supplied from Noord-Holland 


Yes 


3-63 (1-27-10-33) 




No 


Reference 



OR, Odds ratio; CI, confidence interval; QIC, quasi-likelihood under the independence model criterion. 
Number of observations used: 208. Number of levels used: 105 (QIC = 219-1157). 
* See animals at <5 m or touch animals. 



in cattle bulk tank milk was found, confirming wide- 
spread circulation of the bacterium in cattle [36]. To 
further assess the risk for human infection from 
cattle, a similar study addressing the seroprevalence 
and risk factors in dairy cattle farm residents is 
being finalized in The Netherlands. A role for cattle 
in the human infections observed in the current 
sheep farm study, is also supported by the fact that 
the high seroprevalence in sheep farm residents does 
not seem to correspond with the low sheep seropreva- 
lence at the participating farms (<2%). The role of 
specific activities with sheep for the infection risk 
was presumably relatively small, although not absent 
taking into account the significant association between 
human and sheep seroprevalence at the participating 



non-dairy farms. Whether sheep themselves are at 
increased risk for infection because of contact with 
cattle or nearby goat populations is currently 
under investigation. In The Netherlands, a dominant 
C. burnetii genotype was identified in humans, goats, 
and sheep throughout the entire affected area; the gen- 
otype found in cattle appeared to be different [37, 38]. 

Based on the results of the present study, some rec- 
ommendations can be made. First, we want to eluci- 
date the transmission cycle between different species 
of ruminants and farm residents; strains from goat, 
sheep, cattle, and sheep farm residents could be sub- 
typed and compared. Second, more research is needed 
to investigate whether the Blessumer breed is more 
often infected compared to other breeds. Third, 



in this study a high seroprevalence in spouses was 
found (36-3% non-dairy farm spouses, 50-0% dairy 
farm spouses). Therefore, we emphasize the impor- 
tance of the advice that pregnant women should 
avoid contact with sheep during the lambing season, 
and that they should avoid contact with birth products 
of sheep. Currently, the Dutch Health Council is pre- 
paring an advice about vaccination of high-risk pro- 
fessionals, including several farm populations. For 
this advice, they also will take into account the results 
of this study. 

Limitations 

The study of non-dairy sheep farms had a low 
response rate of 9-1%. As reported by several farmers 
not willing to participate, sheep were outside when 
the request to participate was made, and it would be 
too labour-intensive to collect about 60 sheep for 
blood samphng. In addition, this part of the sheep 
industry was not affected by the implemented control 
measures, mainly targeted at farms with dairy sheep 
and dairy goats. Therefore, non-dairy sheep farmers 
might be less motivated to participate compared to 
the small dairy sheep sector, which had a response 
rate of 38%. 

Except for differences in sheep density in the sur- 
roundings and the number of sheep on their farms, 
participating and non-participating non-dairy sheep 
farms appeared to be comparable. As both factors 
were not related to seropositivity, this selective 
response is not thought to be of iniluence on the 
study results, which are therefore considered represen- 
tative for the Dutch professional non-dairy sheep 
sector. 

At 79% of the 119 participating non-dairy farms 
both the farmer and partner participated in the 
study. Therefore, results for the farmers and partners 
are considered representative of the group of farm- 
ers/partners at the participating farms. It was not 
registered how many children aged ^ 12 years lived 
at the participating non-dairy farms, and we cannot 
be absolutely sure that the participating children 
were representative of all children in this age category. 

CONCLUSION 

This study demonstrates that C. burnetii infection is 
common in individuals living and/or working at a 
sheep farm in The Netherlands. Except for their 
sheep, the risk also seems dictated by contact with 
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cattle at present or in the past and by nearby goat 
populations. 
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